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Abstract 

We consider the explanation of both data of the atmospheric neutrino and the solar 
neutrino by the neutrino mass matrix derived from the mixing of neutrinos and gauginos 
in the extended MSSM with an extra U(l) gauge symmetry. This scenario directly relates 
the neutrino mass to supersymmetry. The structure of the mass hierarchy and the mixing 
is determined only by the extra U(l)-charge of neutrinos. Although the model is rather 
simple, it may be able to realize both the small and large mixing angle solutions for 
the solar neutrino problem. In particular, the large mixing angle MSW solution for the 
solar neutrino problem can be easily realized. We discuss the relation between the extra 
U(l)-charge and four solutions for the solar neutrino problem. 



*e-mail:suematsu@hep.s. kanazawa-u.ac.jp 



1 



The atmospheric and solar neutrino observations have strongly indicated the existence 
of the neutrino oscillation and then the non-zero neutrino mass [JE], 0, |3|. This is an only 
known evidence suggesting that there is new physics beyond the standard model (SM). 
Theoretically, the SM has been claimed that there are some unsatisfactory features. Such a 
representative problem is a gauge hierarchy problem. Supersymmetry has been considered 
to be a promising candidate to solve it and its various phenomenological features have 
been studied almost for twenty years Q. However, still now there seems not to be an 
experimental support for it, except for the gauge coupling unification ||. If we can find 
that there are some relations between the supersymmetry and the neutrino mass, it is 
very interesting and might also cause a big impact to the study of the new physics. 

In general the supersymmetric model has a special symmetry called i?-parity defined 
by R p = (— i^ 3B+L+2S , where B is the baryon number, L is the lepton number and S is the 
spin. The ordinary particles in the SM are _R p -even and their superpartners are R p -odd. 
As far as the Z?-parity is conserved, neutrinos cannot mix with neutralinos. However, the 
superpotential can include an i?-parity violating term like e a L a H 2 || [7|, for example, as 
a result of the gravitational effect. Under the existence of this bi-linear f?-parity violating 
term, the vacuum expectation values (VEVs) of sneutrinos is generated and then the 
mixing among neutrinos and neutralinos is induced. The generation of Majorana neutrino 
mass in the minimal supersymmetric SM (MSSM) through this mixing has been discussed 
in a lot of works f| ||, . In most of them only one of three neutrinos is massive at the 
tree level || and then they could be applied only to the atmospheric neutrino problem 



10]. It is necessary to improve this feature in order to explain all of the recent neutrino 



oscillation dataQ. One of the possibilities to make other neutrinos massive is to introduce 
the one-loop effect JT2, 13fl . In such a framework the solar and atmospheric neutrino 
problems have been also discussed [|H|. 

In this paper we would like to propose the generation of the small neutrino mass at 
the tree level based on a typical interaction among neutrinos and gauginos in the usual 
supersymmetric model. The model is very simple and economical. It needs only a small 
extension of the MSSM by an extra U(l) gauge symmetry. There are some works [|14 



on 



1 In this paper we confine our study into the explanation of the atmospheric and solar neutrino data. 
It might be straightforward to extend the model to include a sterile neutrino and make it applicable to 
the LSND data full. 
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the phenomenology including the neutrino mass generation in the supersymmetric model 
with the spontaneous i?-parity violation and the extra gauge symmetry. A new feature of 
our model is that an introduced U(l)-symmetry has generation dependence. The neutrino 
mass matrix produced in this way has noticeable features for the explanation of both the 
atmospheric and solar neutrino data. It might realize the large mixing angle solution for 
them in a natural way if we can assume the spontaneous lepton number violation due to 
the small VEVs of sneutrinos. Our model could be distinguished at least by the existence 
of the additional neutral gauge boson Z' at the TeV region from other neutrino models 



with the large mixing solutions for the solar neutrino problem [15, 16 



In the MSSM there are two neutral gauginos Xw 3 and Ay which are the superpartners 
of the SU(2)i and U(l)y gauge fields. Their interaction with the ordinary left-handed 
neutrinos v a [a = e, fx, r) can be written as M 



Wz v o 



(1) 



If the spontaneous R-parity violation occurs by the VEVs of sneutrinos (D a ) ^ 0, the 
mixing among gauginos and neutrinos is induced. On the other hand, the supersymme- 
try breaking causes the gaugino masses M 2 and M\ for Aw 3 and Ay through a suitable 
mechanism. In the case of M 12 ^> <?i 2(^01)5 we can expect to obtain the small Majorana 
masses of neutrinos as a result of a kind of seesaw mechanism. In fact, if we assume 
(u e ) = (u^) = (u T ) = u, the mass mixing among neutrinos and gauginos can be expressed 



as H> r 



-UATMN ' + h.c.) and 
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where M T = (v a , — ^Avk 3 , iAy). We consider a special case in which every VEV of sneutri- 
nos is equal in order to reduce the number of parameters and simplify the mass matrix of 
the model. Under the assumption of M 12 ^> QxpU we can easily diagonalize the mass ma- 
trix M.. Although there are three light mass eigenvalues, two of them are zeroQ We must 
2 Even if we do not take the VEVs of sneutrinos equal, the situation is the same. In order to improve 
it the one-loop effect has been taken into account in ref . J?], [ijj . 
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extend the MSSM to apply this type of mass generation of neutrinos to the explanation 
of the atmospheric and solar neutrino data. We consider a model which can naturally 
change that aspect at the tree level. 

Let us consider the introduction of an extra U(l)x gauge symmetry specified by the 
following features^. It is non-anomalous and also has flavor diagonal but non-universal 
couplings. The different charge of U(l)x is assigned to the fields belonging to the different 
generationF]. It is assumed to be spontaneously broken at a TeV scale region and then 
we have an additional neutral gauge bozon Z' which is heavier than the ordinary Z . We 
need a new SM singlet chiral superfield whose scalar component causes the spontaneous 
breaking of U(l)x by its VEV. Since a part of the MSSM contents is assumed to have its 
charge in a generation dependent way, it could generate a non-trivial texture in the mass 
matrices. However, we cannot use this symmetry for the Froggatt-Nielsen mechanism [|TJ] 



to induce the hierarchical structure of quark mass matrices since its breaking scale is too 
small as compared to the Planck scale. For simplicity, we do not assign the U(l ^-charge 
to quarks and we consider a model in which only leptons have its charges as 

h ?j, qui), in (-Qi, ~ Qi, ~ Qm)- (3) 

As a result of such a charge assignment u a has a charge q m different from other neutrinos 
vi and vji whose charge is defined by q^. The doublet Higgs superfields are also assumed 
to have no charge so that the ordinary Yukawa couplings of the charged leptons induce 
the texture in the mass matrix such as m IM = m nm = 0. Using this charge, the coupling 
between neutrinos and the U(l)x gaugino is given by iy/2g x J2 a qa (p^XxVa — ^x^aPaj 
as eq. ([I]). We do not consider the kinetic term mixing between the U(l)-gauginos [pDfl . If 
we take this effect into account, off-diagonal elements appear in the gaugino mass matrix. 
In general, the introduction of U(l)x to the MSSM requires addtional chiral superfields 
to cancel the gauge anomaly which causes the non-trivial constraint on the charge as- 
signment. Since its breaking scale is in the TeV region, the non-universal couplings of Z' 
3 The additional U(l)-symmetry is known to appear very often in the heterotic superstring models. It 

can play some useful roles in the supersymmetric models |l7j . 

4 This kind of charge assignment of U(l)x has been discussed in the different scenario to explain the 

small neutrino mass and the proton stability, for example, in |1S| ]. 

5 At this stage we cannot determine to which flavor each v a corresponds so that we used the Roman 

numerals for a. 
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with leptons might impose the constraints on the model through the electroweak preci- 
sion measurements and the flavor changing neutral current (FCNC). We will discuss these 
points later. 

First we focus our attention to the role of U(l)x in the generation of the neutrino mass 
due to the mixing among neutrinos and gauginos. Since the gaugino sector is extended, 
the matrix (^) is modified intof] 

/ x [ a 2 a x b ^ 

M 



m T 
m M 



m 



ci2 ax b 
a>2 a\ b 
\ a 2 a\ c J 



( 



M 



\ 



M 2 







M\ 









M 



(4) 



x J 



where 



V2 



u 



, b = \/2g x q I u and c = s/2g x q m u. We can obtain the light neutrino mass 



matrix from this by using the generalized seesaw formula and it can be written as 



m T M 1 m 



^ m + e 2 m + e 2 m Q + e5 
m + e 2 m Q + e 2 m Q + e5 
i m + e5 m + e5 m + S 2 , 



(5) 



where m , e and S are defined by 



m 



g\u 2 g\u 2 



9xQi u 



9 x (lm u 



(6) 



2M 2 ' 2Mi' v^/.v' VMx' 

The interesting aspect of this mass matrix is that it is defined only by the gaugino mass 
Ma (A = £,X), the gauge couplings g A , the U(l)x-charges q a and the VEV u of sneutri- 
nos. 

We define the mass eigenstates by Mi = (U T ) ia Af a . The diagonalization of the matrix 
(D gives 



U 
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1 


1 




V2 


V2 




cos 8 


cos 8 




V2 


V2 


V 


sin 8 


sin 8 


V2 


V2 



^ 
— sin 9 

cos 9 



(7) 



where one of the mixing angles sin 9 is defined as 

8(m + e5) 2 



sin z 29 



(mo + 2e 2 - 5 2 ) 2 + 8(m + ed) 2 ' 



(8) 



6 In the phenomenological point of view the result obtained in this paper using this mass matrix is 
independent of the mass generation mechanism. If we can find any models giving the neutrino mass 
matrix with the features presented here, the result could be also applied to them. 
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(a, 13) 


(m) 


-AU m U Pl U aj U P] (= A) 




(I, II) 


(1,2) 


cos 2 9 


(A) 




(1,3) 


sin 2 # 


(B) 




(2,3) 


- sin 2 9 cos 2 9 


(C) 


(I, III) 


(2,3) 


2 sin 2 9 cos 2 9 


(D) 


(II, III) 


(2,3) 


2 sin 2 9 cos 2 9 


(E) 



Table 1. The contributions to each neutrino transition process v a — > vp from each sector (i,j) of the 
mass eigenstates. 



The non-zero mass eigenvalues are 

i |(3m + 2e 2 + 5 2 ) =F ^/(m + 2e 2 - 5 2 ) 2 + 8(m + e<5) 2 } . 



^2,3 



(9) 



Here we can consistently assume the charged lepton mass matrix is diagonal^. In this case 
the above mixing matrix U is just the flavor mixing matrix which controls the neutrino 
oscillation. In the following discussion we assume it in the charged lepton sector. There 
is no CP violation in the lepton sector under this assumption. 

It is well-known that the transition probability due to the neutrino oscillation v a — ► z/g 
after the flight length L is written by using the matrix elements of (0) as 

(10) 



f Am 2 - ^ 

V Ua ^ V/3 {L) = 5 af3 - ^Y^UaiUpiUajUpj sin 2 y~^r L 



where Am 2 - 



I 2 

\mf 



In Table 1 we summalize the contribution to each neutrino 
transition mode (a,/3) from each sector (i,j) of the mass eigenstates in our model. As a 
phenomenologically interesting case, we consider the case in which the mass eigenstates 
u 2 and z> 3 are almost degenerate and the mass hierarchy mi m 2 ~ 7713 is satisfied^} If 
we use it to explain the atmospheric and solar neutrino data, the difference of the squared 

7 Although the XJ(l)x constrains the texture of charged leptons mass matrix as mentioned before, 
some additional symmetry might be necessary to make it diagonal. The large mixing between fj, and t is 
disfavored in this model. 

8 This is a well-known reversed hierarchy scenario to the atmosphetic and solar neutrino problems 
However, the absolute value of each mass eigenvalue is smaller than the usual scenario because of 

mi = 0. Every neutrino cannot be a hot dark matter candidate. 
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mass should be taken as 

2 x 1(T 3 eV 2 < Am 2 2 ~ Am 2 3 < 6 x 1(T 3 eV 2 , (11) 

10 -10 e y2 < Am 2 g < 15 x 1Q -4 e y2_ ^ 

The suitable value of Am 2 3 should be chosen within the above range depending on which 
solution is adopted for the solar neutrino problem. We can easily find the simultaneous 
explanation of both deficit of the atmospheric neutrino and the solar neutrino if we identify 
the weak eigenstates (r, /U, e) with (/, II, III). Using this identification, eq. (|7|) is rewritten 
into the MNS mixing matrix in the usual basis as, 






— sin 9 


cos 9 
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cos 9 


sin 9 




V2 


V2 


V2 
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cos 9 


sin 9 
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V2 


V2 


V2 



(13) 



The atmospheric neutrino is found to be explained by — ► v T composed of (A) and (B) 
in Table 1. Here we should note that m\ = and also Am 2 2 ~ Am 2 3 is satisfied. This 
explanation is independent of the value of sin^. On the other hand, the solar neutrino 
is expected to be explained by u e — > (E) and also v e — > v T (D). In both processes the 
amplitude A is | sin 2 29. If sin 2 29 ~ 1(T 2 , it could realize the MSW small mixing angle 
solution (SMA) p2[ . If sin 2 29 ~ 1, it could give the MSW large mixing angle solution 
(LMA), the low mass MSW solution (LOW) and the vacuum oscillation solution (VO) [p^] 
depending on the value of Am^. The CHOOZ experiment constrains a component 
U e i of the MNS mixing matrix in this scenario since the amplitude A of the contribution 
to v e —>■ v x with Am 2 2 or Am 2 3 always contains it. This model is free from this constraint 
since U e \ = 0. 

In Table. 1 an only remaining contribution (C) to — > v T cannot imply any evidence 
in the short baseline experiment even in the case of sin 2 29 ~ 1 since Am| 3 is too small. 
However, this mode might be relevant to the long baseline experiment in the case of 
Am 2 3 ~ 1CT 4 eV 2 which corresponds to the LMA solution of the solar neutrino deficit. 
We show the effect of the mode (C) on the Piv^ — > v x ) in Fig. 1. The dashed line comes 
from the modes (A) and (B). This corresponds to the ordinary two flavor oscillation. The 
thick solid line is the one which is obtained by taking account of the contribution of (C). 
In the thin solid line the contribution of (D) which corresponds to v x = v e is also taken 
into account. This shows that it might be possible to discriminate this model from others 
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Fig. 1 The transition probability Piy^ — > v x (^^)) as a function of the flight length L. We assume 
E =1 GeV, Am 2 3 = 3.5 x 1(T 3 eV 2 and Am| 3 = 1(T 4 eV 2 

in the long baseline experiment such as L > 2000 km. Moreover, this model may be 
expected to have another experimental signature in the neutrinoless double /3-decay . 
Using eq. (|i~3D, the effective mass parameter which appears in a formula of the rate of 
neutrinoless double /3-decay is estimated as 



(m,2 sin 2 6 + cos 2 6^ ~ m^. (14) 



Thus \m ee \ takes the value in the range of 0.04 - 0.08 eV which is independent of the value 
of sin 9, that is, the solution of the solar neutrino problem. The value seems to be within 
the reach in the near future experiment. 

In order to see whether the values of the oscillation parameters can be successfully 
realized, we need to study them by using the parameters Ma, g A , q a and u. In the usual 
soft supersymmetry breaking scenario the gaugino mass is universally produced as Mo 
at the unification scale. Its low energy value is determined by the renormalization group 
equations (RGEs). If we use the one-loop RGEs, the gaugino mass at a scale \i can be 
expressed as 

MM = ^fl£(/*), MM = l^g'M, (15) 
9u 3 9u 

where we assume the gauge coupling unification at the scale My and define a value of the 

gauge coupling at M\j as g u . It is not unnatural to consider the gauge coupling of U(l)x 

and its gaugino mass to be the same as the ones of U(l)y f\. If we adopt this simplified 

9 It might be satisfied if U(l)x is unified into a simple group together with other SM gauge group like 
the SO (10) and E6 models. Also in the superstring context the freedom of an Abelian Kac-Moody level 
could make it possible. 
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2 

possibility, we can find that m 2j 3 and sin 2 26 can be written by using only q a and j^u 2 as 



™ 2 ,3 = ~ I 2 + 2 ?/ 2 + =F y (2 + 2^ + )2 _ ^ _ gJa j 2jX (16) 

sin 2 2^ = ? 2 8(2 2 + 3 ^^ )2 ^. (17) 

(2 + 6g 2 -3g 2 ) 2 + 8(2 + 3 g/ ^) 2 1 ' 

The structure of the mass spectrum and the flavor mixing is controled by the U(l)x- 
charge. The gaugino mass Mq and the VEV u of sneutrinos are relevant to the mass 

2 

eigenvalues only in the form of an overall factor j^u 2 . In order to realize the value of (|TT| ) 

2 

which is required by the atmospheric neutrino deficit, j^-u 2 needs to be in the range of 
0.017 eV < |U 2 < 0.023 eV Q If we take M ~ 100 GeV and g v -0.72, for example, it 
shows that the sneutrino VEV should be u ~ 60 - 70 keV. The remaining freedom which 
we can use to explain the solar neutrino deficit is the U(l)x-charge of neutrinos. In Fig. 
2 we plot the value of the U(l)x-charge and the corresponding oscillation parameters 
which are suitable for the explanation of the solar neutrino deficit. The figure shows that 
the reasonable value of the U(l)x-charge can realize all the well-known solutions for the 
solar neutrino problem. As is easily seen from eq. (16), Am^ is propotional to the part 
of square root in the expression of m 2j 3. It means that the solution with the smaller 
Am 2 ^ requires the finely tuned U(l)x-charge^|. As a result, the SMA, the LOW and 
the VO seem to require the finer tuning of the U(l)x-charge than the LMA. It might be 
favorable for this model if we take seriously the recent Super-Kamiokande analysis of the 
solar neutrino ||. Anyway it is interesting that every solution is obtainable within this 
simplified framework. 

Now we comment on some remaining important points. The first point is whether the 
anomaly-free U(l)x can realize the solutions in Fig. 2(a). We show it by using a simple 
example presented in (|3|). In order to cancel the gauge anomaly we introduce additional 
chiral superfields with the non-zero U(l)x-charge: 2o, l±i and 4(lo), where we indicate 
the representation under SU(2)i,xU(l)y. The U(l)x-symmetry is assumed to be broken 
by the VEVs of Iq's. The 2 and 1±\ should be massive at the TeV scale through the 
VEVs. If we impose this requirement with the anomaly-free conditions on the U(l)x- 

10 In the estimation we have already taken account of the effect coming from the U(l)js:-charge depen- 
dence. However, its effect is not large and only a factor of 0(1). 

ll lt may be useful to note that the complete degeneracy 7712 — 7713 occurs in the case of q r — l/y/3 and 

q B = -2/V3, which satisfy Y< a Qa = 0. 
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Fig. 2 (a) Scatter plots of the U(l)js:-charge of neutrinos which can give the solution to the atmospheric 
and solar neutrino problems. The dot-dashed, solid and dashed lines are the ones corresponding to 
Q = 1.04, 1.05 and 1.06, respectively, (b) Oscillation parameters for the solar neutrino which are realized 
by the U(l)x-charge of neutrinos shown in (a). 



charges, we find that the charges of all fields can be expressed by using the charge Q of 
one of lo's. In Fig. 2(a) we plot the values of (qi, qm) of such solutions by taking suitable 
values of Q. From it we find that our scenario could have the solutions consistent with 
the above mentioned conditions. Although we take a very simple example to show it, it 
would be possible to construct more elaborate examples. 

The next point is the phenomenological constraint on the U(l)x since we have a Z' 
at a scale not far from the weak scale. The constraint on the Z' comes from the direct 



search [£5|] and the electroweak precision measurements [^g, [I7| . In general these seem 
to require the mass bound vnz> > 600 GeV, which can be consistent with our assumption 
of the XJ(l)x breaking scale. Another constraint could also come from the FCNC in the 
lepton sector such as the coherent fi-e conversion, r — > 3e, 3/i and [i — > etc., since 
we assume the flavor diagonal but non-universal couplings of Z' to leptons. The detailed 



analysis for such a issue has been done in |27]] and we can use the discussion there. In our 
model the Z 1 interaction term can be written in the mass eigenstates as, 



£z> = -gi 



9x 
9i 



^cosCJ ( x) - sin ^fl) 



Z' 



10 



Bg = V^diag(gar,g f ,«r)^, (18) 

where £ is a Z-Z' mixing angle and is the unitary matrix to diagonalize the mass 
matrix of ip. The FCNC in the lepton sector causes the strong constraint on B^' R as it 



is discussed in P7[ . Since the U(l)x-charge of fi and r is equal in our model, B 2 ^' R = 0. 
Moreover, if we remind that m eA1 = m er = is satisfied in the charged lepton mass matrix 
because of the constraint of the U(l)x charge, we find that b\ 2 r = B^' R = 0. Then 
the non-universal couplings of Z' do not induce the observable effect on the FCNC. The 
addtional contribution to the FCNC could be caused by the extended neutralino sector. 



The study in that case has been done in [20] for the models which can be embedded in E 6 . 
The effect is found not to be so large as far as we take parameters in such a way that the 
MSSM satisfies the experimental bound. Although the U(l)x-symmetry is different from 
this case, the qualitative feature is similar also in the present model. Thus as far as we 
consider a TeV scale Z' with the small mixing £ < 1CT 3 , the constraint from the precision 
measurement and the FCNC of the lepton sector seems not severe to the present model. 

Another important point is the origin of the small VEV u of sneutrinos. As mentioned 
in the previous part, it should be around O(10 2 ) keV which is much smaller than the weak 
scale. In the MSSM there are arguements on the lepton number violation due to the VEVs 
of sneutrinos in the vicinity of the weak scale p8[ and also some authors point out that 



the neutrino mass produced by them can be sufficiently small |13| , Q. However, in our 



scenario we need much smaller VEVs of sneutrinos than the weak scale. As one of the 
possibilities, we might consider that such small VEVs of sneutrinos could be obtained if 
there is bi-linear i?-parity violating terms eL a H 2 with the sufficiently small e 0. We can 
check it briefly by minimizing the scalar potential under the assumption that (^1,2) can be 
treated as constants. In this case the value of u derived from the potential minimization 
condition can be approximately written as 

u 6e^( Hl ) + B e (H 2 )) 

3(3i 2 + <? 2 2 )«#i> 2 - (#2> 2 ) + 2g 2 x (Ea q^iqiW + g 2 (F 2 ) 2 ) + 6m 2 ' 1 ) 

where B 6 is a soft supersymmetry breaking parameter related to the eL a H 2 terms and 
12 We assume = e. The smallness of e also causes a new hierarchy problem. To resolve this problem 
we might be able to use the similar mechanism to the solutions for the [i problem (29), although the 
relevant energy scale needs to be scaled down by some orders of magnitude. 
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m 2 is the soft scalar mass of sneutrinos, which is assumed to be universal. From this 
we find that the sufficiently small u can be obtained as far as e is small enough and 
the /z-parameter, B t and m 2 take the values of the order of weak scale. However, the 
eL a H 2 terms generate the y a -H 2 mixing which affects our analysis. In order to escape the 
influence due to the e terms, ^ might be required to be sufficiently larger than -. We 
need to check whether these conditions are satisfied at the true vacuum with the equal 
VEVs of sneutrinos taking account of the radiative correction. Including the study of 
other possibility to realize the small u, we would like to leave it in the future publication 
instead of going further into it here. 

In summary we have proposed the scenario for the mass and the mixing of neutrinos in 
the super symmetric model with an extra U(l)x-symmetry. The scenario is based on the 
mixing among neutrinos and gauginos. In this model we could obtain the non-zero mass 
eigenvalues for neutrinos at the tree level. The atmospheric and solar neutrino deficits 
can be simultaneously explained by the reversed mass hierarchy scenario. In particular, 
every known solution for the solar neutrino problem can be realized only by tuning the 
U(l)x-charge of neutrinos. It is interesting that the large mixing angle MSW solution 
can be easily realized as compared to other solutions. The neutrinoless double /3-decay 
might be accessible if the experimental bound is improved to the level of \m ee \ ~ 0.04 - 
0.08 eV. There remain some unsolved problems related to the VEVs of sneutrinos and 
the supersymmetry breaking parameters. Further investigation of these problems seems 
to be necessary to see whether our model works well in a realistic way. Finally we would 
like to stress that the texture of the present mass matrix might be scaled up into the 
high energy region. In that case it may be regarded as the usual seesaw model with a 
special Dirac mass matrix. Various conditions on the U(l)x-charge etc. in the present 
scenario could be replaced into the ones for Yukawa couplings and so on. The study in 
this direction also seems to be interesting. 

I would like to thank M. Tanimoto for useful comments related to the CHOOZ constraint. 
This work is supported in part by the Grant-in-Aid for Scientific Research from the 
Ministry of Education, Science and Culture (No. 11640267). 
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